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Abstract. I discuss many-body models for interacting fermions in two space 
dimensions which can be solved exactly using group theory. The simplest example 
is a model of a quantum Hall system: 2D fermions in a constant magnetic field and 
a particular non-local 4-point interaction. It is exactly solvable due to a dynamical 
symmetry corresponding to the Lie algebra gl^^ ® gloo- There is an algorithm to 
construct all energy eigenvalues and eigenfunctions of this model. The latter are, in 
general, many-body states with spatial correlations. The model also has a non-trivial 
zero temperature phase diagram. I point out that this QH model can be obtained from 
a more realistic one using a truncation procedure generalizing a similar one leading to 
mean field theory. Applying this truncation procedure to other 2D fermion models I 
obtain various simplified models of increasing complexity which generalize mean field 
theory by taking into account non-trivial correlations but nevertheless are treatable by 
exact methods. 



1. Introduction 

In this paper I present a class of exactly solvable many-body models which describe 
two dimensional (2D) interacting fermion systems with spatial correlations. By 'exactly 
solvable model' I mean that there is an algorithm to construct all eigenvalues and 
eigenstates of the many-body Hamiltonian explicitly. By 'correlated' I mean that the 
eigenstates, and in particular the groundstate, are not Slater states in general. I 
introduce these models by discussing a simple example which I call QH model. This 
is a model of 2D electrons in a magnetic field and interacting with a particular 4-point 
interaction p. I then propose a novel interpretation of this model allowing for natural 
generalizations and leading us to a large class of exactly solvable models of interacting 
fermions. As I will explain, all these models can be obtained from more realistic ones by 
a truncation procedure which generalizes a similar one leading to mean-field theory. As a 
motivation, I will discuss two such models which played an important role in condensed 
matter physics: the BCS model which solved the problem of superconductivity, and the 
Hartree-Fock model providing a basis for understanding magnetism in metals. These 
models can be obtained from a more general one by discarding all but particular kinds 
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of terms in the interactions. The Hartree-Fock model is obtained by keeping only the so- 
called Hartree- and Fock terms, i.e. the terms describing two-body scattering processes 
where the fermion momenta are left unchanged and are exchanged, respectively; see the 
Fig. in ()36|) below. In 2D there are also mixed scattering terms which are Hartree-like 
in the x- and Fock- like in the ^/-component of the momenta and vice versa; see Fig. fl54|l . 
I observe that the interaction in the QH model is a sum of such mixed terms, and I 
propose generalized models where, in addition to Hartree- and Fock terms, also mixed 
terms are included. These models are interesting since they include non-trivial spatial 
correlations but still can be treated by exact methods. 

I only mention in passing some physics motivation (and refer readers interested in 
more details to a review [2'.). There are several fascinating phenomena which have been 
discovered in complex transition metal oxides in recent years. A prominent example is 
high temperature superconductivity. 2D models of Hubbard-type are generally beheved 
to account for these phenomena, but these models have been found to be very difficult: 
despite of much work they are still poorly understood. I therefore believe that the 
simplified models proposed here will be useful in this context. 

It is important to note that the exact solution of Hartree-Fock type models is 
equivalent to a mean field approximation of a more complicated model (I will discuss 
this in more detail in Section [3.5.2|) . and the same is true for the BCS model. § We 
thus propose that our 2D models provide a generalized mean field theory taking into 
account non-trivial spatial correlations. The latter is somewhat more complicated but 
still can be computed by exact methods. To be more specific: Standard mean field 
theory corresponds to minimizing a functional depending on a small number of variation 
parameters, whereas the generalized mean field theory is given by a matrix model. E.g. 
the grand canonical partition function of the QH model can be written as a hermitean 
matrix integral 



with a function F only depending on the eigenvalues of the L x L hermitean matrix X, 
and the thermodynamic limit corresponds to L ^ cxo. Much is known how to treat such 
matrix models exactly (see e.g. [3] for an introduction to this field). I believe that it is 
possible to find exact expressions for these partition functions. 

The simplified models discussed in this paper are such that eigenstates and 
eigenvalues can be found exactly using group theory. Models of this kind are standard 
in nuclear physics |1] but have been used surprisingly little in condensed matter physics. 
An exact solution of a BCS model j3] found by Richardson a long time ago [H] went 
unnoticed in condensed matter physics until recently [7 . This model and some recent 
generalizations fS] are similar to the models proposed here in that all of them can be 
solved using group theory. In fact, the model solved in ^ is similar to what we call 

§ I believe this is the reason why these models are no longer widely known: there are other, more 
efficient methods to derive mean field theory. 




(1) 
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Hartree-models in that the corresponding Lie groups are finite dimensional (SU(2), e.g.), 
whereas our 'mixed' models are related to infinite dimensional Lie groups. 

The plan of the rest of this paper is as follows. In the next section we discuss 
the QH model. I point out various remarkable properties of this model and present its 
solution obtained in Section 3 contains a discussion of general interacting fermion 
systems and the truncation procedure leading to simplified models which can be solved 
exactly. The 2D Hubbard model is a special case where our truncation procedure leads 
to particular simple models. We end with a short summary in Section 4. 

An outline of this article appeared in jH]. 

2. The Quantum Hall model 

In this Section I discuss a model of interacting fermions proposed and solved in p]. 

2.1. Definition and physical interpretation 

We start with the Landau Hamiltonian 

hs = (-19. + By)'' + i-idy - Bxf (2) 

describing the motion of a charged particle in the (x, ?/)-plane with a constant magnetic 
field of magnitude -B > perpendicular to this plane. We choose units such that the 
mass and charge of the electron are 1/2 and 2, respectively. A microscopic model for a 
quantum Hall system|| is given by the many-body Hamiltonian H = Hq + V where 

Ho= f d2xV'^(x)(/iB-/i)^(x), (3) 

with fermion creation- and annihilation operators ip and tp"^ obeying the usual canonical 
anticommutator relations, {?/'(x), (x')} = 5^(x — x') etc.; x = {x,y), and the real 
parameter is a chemical potential. The interactions V is of the following two-body 
type, 

V = j d2xi...d2x4t;(Xi,...,X4) : V^t(Xi)V^(x2)V^t(x3)^(x,) ; (4) 

with the dots indicating normal ordering as usual. ^ The interaction vertices v of main 
physical interest for quantum Hall physics are as follows, 

f (Xi, . . . , X4) = 5^(xi - X2)5^(X3 - X4)iy(xi - X3) (5) 

with W the two-body potential. For example, Wi^x) = const./ |x| corresponds to a 3D 
Coulomb repulsion, but this choice leads to a model which cannot be analyzed without 
approximations. One might hope to obtain a simpler model by taking a fully local 

II I do not write 'quantum Hall effect model' since, firstly, such a model should also include impurities 
and, secondly, I do not know if the interactions I propose below can account for the QHE. Note also 
that I ignore spin. 

1 : V'ni)V'(2)^U3)^(l) : = (3)^1(1)^(2)^(4) 
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interaction iy(x) = 5'5^(x — y), but this model is trivial since the fully local interaction 
vanishes due to the Pauli principle, V = 0. 

The exactly solvable QH model can be obtained as a deformation of a fully local 
interaction which is such that it becomes non-trivial. In Fourier space this deformed 
interaction is given by 

K3H=-^ / d'kl...d2k4<52(ki-k2 + k3-k4) X 

^g-ie(kiAk2+k3Ak4) . ^t(ki)^(k2)V^^(k3)V^(k4) : (6) 

with the hat indicating Fourier transform, k A q = kxQy — kyQx, and 6 the deformation 
parameters (i.e. ^ = corresponds to the fully local interaction). At first sight this 
interaction might appear somewhat unusual, but it has several remarkable mathematical 
properties which I will discuss below. These mathematical properties were my original 
motivation to consider this model pp. In the next Section I will discuss other interacting 
fermion models which suggest an alternative physical interpretation. This discussion will 
lead us to various interesting generalizations of this QH model which, as I believe, are 
promising candidates for understanding 2D correlated fermion in a precise mathematical 
framework. 



2.2. Mathematical properties 

The QH model defined above has various interesting mathematical properties. Firstly, 
the interaction in Q can be obtained from the fully local interaction by replacing the 
pointwise product of fields by the so-called Groenewold-Moyal star product, 

VQn = 9o[ d^x :7/>t(x)^V>(x)*^t(x)^^(x) : (7) 

with -k the associative product defined such that x-ky = xy — i9, y-kx = xy + i9, x-kx = x^ 
and y-ky = y"^. Indeed, from these rules one can deduce that the T<r-product of two plane 
waves is e^'^^'^^e^'^''' = e^^'^+'^^^'e"^^''^'', and with that one can compute Eq. ((7j) by inserting 
the Fourier transforms of the fields and obtain (jH}. As can be seen by searching the 
arXive, field theories with such T^r-interactions received much attention in the particle 
physics literature theory recently (they are usually referred to as 'noncommutative field 
theories'). Secondly, the interaction Vqh has the remarkable property that it looks the 
same in Fourier- and position space: Computing the inverse Fourier transformation one 
finds that Vqh can be written as in Eq. with the interaction vertex ^01 

t;(xi, . . . , X4) = - X2 + X3 - ^^)e--^(-i^--+-^^-^y<^ . (8) 

Since the Landau Hamiltonian obviously looks the same in position- and Fourier space 
as well, the Hamiltonian H = Hq + Vqh is invariant under the following duality 
transformation [TU] , 

z^(x) V'(x) = Bi:{B^) 
B ^ B = B 
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e -^6 = -1/(63^) 

90 ^~9o = 9o/\Be\. (9) 

Note that this duahty transformation involves Fourier transformation, ^p{'x) —>■ 'ip(k), 
followed by a re-scaling of the Fourier variable to correct for the length dimensions, 
k — > k/i? = X (since B has the dimension length'"^). Thirdly, as I will discuss in 
more detail below, this model is exactly solvable for 6 = ±B^^ [T\. It is interesting to 
note that these exactly solvable cases are mapped onto each other under the duality 
transformation. 

We can actually give the solution of a more general model where /ig in Eq- & is 
replaced by {1 — a)hB + ah_B with < a < 1. This generalization means that we allow 
for a confining electric potential oc (x^ + y^) in addition to a magnetic field. 



2.3. Integrability 

In the following we assume 6 = 1/B. 

As discussed in most quantum mechanics textbooks, the eigenf unctions 
of the Landau Hamiltonian in Eq. Q are labeled by two positive integers £, m with 
eigenvalues depending only on one of these quantum numbers, hB(p£m = ^£(pim where 
€£ = 4:B{i—^) (this is the well-known degeneracy of the Landau problem). The quantum 
numbers i and m can be interchanged either by complex conjugation, 0£m(x)* = ^^^(x), 
or by changing the direction of the magnetic field, h_B(pem = ^m4>tm- Note that changing 
the sign of the magnetic field is equivalent to a parity transformation, x — x. 

These eigenf unctions provide a complete orthonormal basis in the 1-particle Hilbert 
space L^(R^), and thus the many-body Hamiltonian Hq in Eq. Q can be diagonalized 
by expanding the fermion field operators in this Landau basis, '?/'(x) = J2em'^^"i4>im{^) 
and similarly for V'^^)- This yields Hq = J^emi'^i ~ l^)^\m^tm- Since the (p^rn are 
common eigenfunctions oihB and /i__b, we can actually generalize the model by replacing 
/i^— >(1 — a)/iB + ah^B and still obtain a diagonal many-body Hamiltonian, 

with = {\ — a)ei and = acm — fJ'- The fermion field operators cj^ and c^^ obey 
the canonical commutation relations {c^,i^,cl,^,} = Se^e'Sm,m' etc. The model is defined 
on a Fermion Fock space with a vacuum Q such that CimQ = for all i, m. 

The crucial fact making the QH model 'special' are the following remarkable -k- 
product relations of the Landau eigenfunctions [Ij, 

1 

(pim'^ <Pi'm' = \ -7— ^m,&4>tm' ^OI 9 = —. (11) 

V 47r B 
Using that and other properties of the Landau eigenfunctions, a straightforward 
computation brings the interaction in ((7j) to the following simple form. 
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where g = goB/An. 

It is interesting to note that the QH Hamiltonian can be written in terms of the 
following fermion bilinears, 

e m 

This allows us to write 

Hq — ^2 EmPmm + Empmm (14) 
m 

and 

Vqh = 5' XI • Pmf-P^'m ~f XI • P(-'^P^(- '■ ■ (1^) 

t,m Tn,£ 

If we regard the operators p£rn as elements of an infinite matrix p and similarly for p we 
can write this in the following suggestive matrix notation, 

Ho = Tt{Ep + Ep), V(^n = g: Tr(p2) := -g Tr{p^) ; (16) 

with E = diag(i?i, E2 . . .) and similarly for E. 

Using the fermion anticommutator relations we derive the following commutator 
relations 

\Pim,Pem'] — ^m,t'Ptm' — ^m! ,iPt'm (17) 

and similarly for the p's, and [permPem'] — 0. Moreover, 

pL = Pme (18) 
and similarly for the p's. We thus see that the operators p and p represent the Lie 
algebra ® gl^^. It is not difficult to check that the interaction Vqh represents a 
Casimir element in this Lie algebra (i.e. it commutes with 'everything') and the free 
Hamiltonian Hq represents an element in the Cartan subalgebra (i.e. is a sum of mutually 
commuting terms). However, the Lie algebra gl^^ ©gloo does not commute with H: one 
can check the following commutator relations, [H,pui] = {E^ — E^i)pu' ■, and similarly for 
p and E. Thus this Lie algebra defines a dynamical symmetry of the QH model. This 
suggests that the model is integrable. We will show this by presenting an algorithm to 
construct all eigenstates and eigenvalues of H below. 

2.4- Solution 

A complete orthonormal basis in the Fock space of the model is given by the states 

with = 0,1,2,... the fermion number. We find it convenient to introduce a 
regularization by restricting the quantum numbers to £j, ruj = 1,2 ... ,L with a finite 
cut-off L. This has a natural physical interpretation: we add a box-like potential 
confining the QH system to the disc + < const. L^/B. This reduces the dynamical 
symmetry to gl^ffiglL, and the fermion Fock space becomes finite dimensional, dim= 2^^ . 
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Straightforward computations using the canonical anticommutator relations show 
that the interaction acts on these states as follows, 

: Tr(p2) : |iV) = 2 T^,,^\N) ^ C^\N) (20) 

with T(jj) transpositions interchanging rrii and mj, 

Tiij)\N) = cl^^ . . . c\^^ . . . c\^^ ■ ■ ■ c\^^^n, i < j. (21) 

Note that T defines a representation of the permutation group Sn, and Cn = J2i<j '^{ij) 
commutes with all permutations Tp. All |A^) are, of course, eigenstates of the free 
part of the QH Hamiltonian, Hq\N) = So\N) with Sq = Yl!j=i{^i-3 + ^mj)-, and due 
to our special form of £q this free part commutes with all permutations Tp. Thus the 
eigenstates of H can be constructed by solving the following eigenvalue equation, 

Cn apTp = v apTp. (22) 

Each solution of i2^) provides and eigenstate ip and corresponding eigenvalue S of the 
QH Hamiltonian H as follows, 

N 

^ = J2apTp\N) and £ = ^{E,^ + E^^^) + 2gv. (23) 
P i=i 
Eq. can be solved by constructing and diagonalizing A^! x A^! matrices representing 
Cn- To appreciate the difficulty of this problem it is instructive to do this for A^ = 2, 3. 
Obviously such brute-force approach is restricted to rather small A^. Fortunately for us, 
Eq. (j22|) corresponds to a classical problem in group theory which was solved a long time 
ago (the group theory results used in the following are discussed in most advanced group 
theory books; for more details and specific references see |I]): The representation T is 
(essentially) the regular representation, and solving Eq. amounts to decomposing 
T into irreps. The irreps of are well-known. They are one-to-one correspondence 
with partitions [A] = [Ai, A2, . . . , Al] of A^ where the A/s are integers summing up to A^ 
and obeying Ai>A2>--->Ai^>l. The possible eigenvalues v in Eq. are given 
by the value of the class operator Cn in the irreps [A], 

L ^ L 

v[x^ = Y.-\^{K + l) - Y^iK (24) 

i=l i=l 

and its multiplicity is the square of the dimension the irreps [A]. E.g. for A^ = 4 the 
eigenvalues c of C4 are 6(1), 2(9), 0(4), —2(9), —6(1) with the numbers in the parenthesis 
giving the multiplicities. The corresponding ^pOpTp is equal to the Young operator 
3^[A] associated with [A]. This gives the following exphcit recipe for constructing the 
eigenstates associated with [A]: For sets of K < N distinct integers {11,12, ... .ir) , 
I < ij < N , define (anti-)symmetrizers -S'^^^^^^ = J^peSK^^^^^^^'^p'^^ '^'^ acting 
on states \N) as described above, with the integers ij indicating the positions of the 
rui. on which the permutations act. Write the integers 1,2, ... ,N into the boxes of the 
Young tableaux associated with [X], take for each row one symmetrizer with indices 
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given by the numbers in the boxes of the row, and similarly for each column take an 
antisymmetrizer. Then 

^ = 3^[A]|iV) = n'^-wsH'^c-olumnslA^)- (25) 

For example, the 10-particle eigenstates associated with [A] = [5, 3, 2] are 

= '5lt2,3,4,5'^(t7,8'^!tl0'^l,6,9'^2,7,10'^3,8l^ = ^0)' (26) 

and £ = J2]=ii^^j +Ek^) + 14:g are the corresponding eigenvalues. This gives all possible 
eigenvalues and corresponding eigenstates of the QH model. It is worth mentioning that 
there is another representation T of the permutation group Sn acting on the indices ij 
of the states |A^) in (|Tn|l . The two representations are, however, related: interchanging 
rrii rrij and £{ ^ ij together is the same as interchanging two fermions and only gives 
a minus sign, and thus TpTp = (— l)'-'^' for all P E Sn- 

It is interesting to note that the model remains solvable if we add to H the following 
interaction terms, 

Vr = VlmpUPmm + V£mPuPmm + WemPuPmm (27) 

with V, V and w arbitrary parameters: the eigenstates ip in Eq. ()23|) remain the same, 
and the eigenvalues £ are changed by adding ^(fmimj + ve-e^ + Wrmij)- We will come 
back to this in Section 13.31 

2. 5. Physical properties 

We found all energy eigenstates and eigenstates of the QH model. However, this is only 
a first step to understand the physics described by this model: much work remains to 
be done. In this Section I describe a few results and mention interesting open problems. 

2.5.1. Correlations. The |A^) in Eq. (fTTH) are Slater states: they are eigenstates of 
the free Hamiltonian in Eq. (fTUI) and describe fermions without any correlations. In 
the next Section we will discuss the so-called Hartree-like models. These are exactly 
solvable models with interactions but such that the eigenstates are still Slater states: the 
interaction only modifies the energy eigenvalues but the eigenstates remain uncorrelated. 
The QH model solved above is remarkable in that its eigenstates are particular linear 
combinations of Slater states. These eigenstates tp still are characterized by pairs 
of quantum numbers {ij,mj), but they are complicated linear combinations of Slater 
states |A^) where the ij^s and rn^'s are distributed over the fermions in many different 
ways. The particular distributions which yield eigenstates are labeled by partitions [A], 
and the energy of the state depends not only on the quantum numbers {ij, rrij) but also 
on [A]. It is impossible, in general, to write these eigenstates ip as Slater states. It 
would be very interesting to compute 4-points Greens functions and thus understand 
the correlations quantitatively. 
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2.5.2. Ground state and partition function. A first step in understanding the physical 
properties of the QH model is to determine the ground state at fixed fermions density 
u = N/L and in the thermodynamic hmit L oo. For that one needs to find the 
eigenstate ip at fixed, large particle numbers N such that the corresponding eigenvalue 
S is minimal. This is a non-trivial problem whose solution I only know in one simple, 
special case. 

The main difficulty in determining the groundstate of the QH model from the 
solution given above is that we actually found too many eigenstates: For each set of 
quantum numbers {ij,mj) we constructed A^! eigenstates, but these can be linearly 
independent only if all £j's and all mj's are distinct. For example, if all m/s in |A^) 
Eq. p9|) are the same, then |A^) is only non-zero if all f^'s are distinct (Pauli principle). 
Moreover, in this case all permutations act trivially, Tp\N) = \N), and thus there is 
only one non-zero eigenstate, namely ip = \N), corresponding to [A] = [A^] (the Young 
tableaux with only one row). The opposite extreme is when all f^'s are the same. In 
this case all m^'s need to be distinct, and the only non-zero eigenstate corresponds to 
[A] = [1^] (the Young tableaux with only one column). In general, degeneracies of 
the ij's and rn^'s will reduce the multiplicities of the eigenstates given above. Many of 
these multiplicities are zero due to the Pauli principle, and to find the minimum energy 
eigenvalue with non-zero multiplicity can be difficult. This seemingly technical point 
has important physical implications: it can lead to a frustration in the system, and thus 
the groundstate can change drastically as the model parameters are varied. Below I will 
illustrate this by discussing the ground state of the QH model as a function of g/B. 

More generally one would like to compute the grand canonical partition function 
Z = Tr exp(— /3iJ). The computation of Z seems like a highly non-trivial problem, but 
I believe it has a beautiful solution, at least for the case where E = 0: as already 
mentioned, one can represent 2 as a matrix model with an external field E, and 
even though the resulting matrix model does not seem to be a 'standard' one it seems 
computeable (it generalizes a matrix model recently solved in fl]). I now sketch how 
to represent ^ as a matrix model. We start with a standard representation of Z as 
a fermion functional integral, i.e. an integral over r-dependent Grassmann numbers 
where < r < /3 is the Matsubara time; see e.g. p2j. One then introduces a Hubbard- 
Stratonivich field Y{t), < r < /?, which allows to integrate out the fermions. The 
fermion integral yields a functional determinant det(i?®/ + /® [9t- — iF (r)]). The field Y 
are hermitean Lx L matrices and periodic, Y{0) = Y{j3). Thus Z becomes a hermitean 
matrix path integral. One then can change variables to X = U (t)^^ [1(9,-+!^ (r)]f/(r) with 
f/(r) unitary matrices and periodic and X a hermitean constant matrix (independent 
of r). The integral over U{t) can be done, and thus one can represent Z as an integral 
over the hermitean matrices X in an external field E. 1 believe it is possible to evaluate 
this matrix integral in the large-L limit by generalizing the computation of a similar 
integral in fTI] . 
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2.5.3. Zero temperature phase diagram of the QH model. I now discuss the ground 
state of the QH model 

H = J2 ^^(^ - 1) + ^ I] : P^"^Prni ■ (28) 

m>l £,m 

as a function of the couphng parameter g/B [B > 0, g arbitrary real, and we set 
= 1/2 for convenience). As discussed, the energy eigenvalues are sums of two terms: 
the kinetic energy £q = J2j'=i^^i^j ~ 1); energy Scott = '^gv[x] due to the 

interaction. We consider the case N < L, i.e. the fermion filling factor u = N/L is 
less than 1. Then the kinetic energy is minimal if all fermions have the same quantum 
number rrij = 1. In this case |A^) in Eq. ()19|) is also an eigenstate of the interaction 
Vqh, and the corresponding eigenvalue is Scott = gN{N — 1). This obviously leads to a 
minimal total energy S^^] = gN{N — 1) if g < 0. However, for positive g it is preferable 
have some fermions with rrij > 1: this increases the kinetic energy but allows to decrease 
the correlation energy. In particular, for very large positive values of g the correlation 
energy will dominate. The minimum value Scon = —gN{N — l) of the correlation energy 
is for [A] = [1^], but we need all m/s to be different to get a non-zero such eigenstate. 
The minimum kinetic energy possible then is for rrij = j. This gives the total energy 
£^[ijv] = {2B — g)N{N — 1), which is lower than E^^t^ for g > B. One can show that the 
latter states are groundstate for g > 2B. The groundstates in the immediate regime 
< g < 2B can be found by solving a variation problem. One finds that they are given 
by partitions 

Aj = non-negative integer close to D — Ki, < i < D/ K (29) 

with K = {2B — g)/g > and D = v2NK (the corresponding Young tableaux 
approximates a right triangle with area N and slope K). One thus finds that the 
QH model has a nontrivial zero temperature phase diagram with zero temperature 
phase transitions at g = and g = 2B. The groundstates are Slater states for g < 
and g > 2B but are correlated states for < g < 2B. All groundstates are highly 
degenerated for 1 <^ <^ A, and the degeneracy (entropy) varies significantly in the 
region < g < 2B. 

3. Interacting fermion systems 

In this Section I will discuss several well-known exactly solvable models of interacting 
fermions which play an important role in condensed matter physics. This discussion 
suggests an interesting physical interpretation of the QH model discussed in the previous 
section. This interpretation will lead us to propose a general strategy to obtain 
'interesting' simplified models for 2D correlated fermions. 

The general class of interacting fermion models which we consider describe fermions 
with 1-particle states labeled by a large but finite number of different quasi-momenta k 
and a spin index a assuming at most two different values ('at most' since we also discuss 
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models for spinless fermions). The models are defined by an Hamiltonian H = Hq + V 
with the free part 

^o = E(^k-/^)4A„ (30) 

where is the dispersion relation and /i the chemical potential. We assume that the 
interaction is of the following form, 

^ = X] X^%i...k4 • '^ki,a'^k2,a'^k3,/3'^k4,/3 • (31) 

a, 13 

where V]^_^ is the interaction potential in Fourier space. The cj^"^^ are fermion operators 
obeying the usual anticommutator relations 

{4,a'%,/3} = -^k.k'^a./S (32) 

etc. The model is defined on the Fermion Fock space with a vacuum Q such that 
Ck,o^ = for all a and k. The primed sum is over all ki, . . . , k4 restricted by momentum 
conservation. In simple cases the latter amounts to ki — k2 + ks — k4 = 0, but if the 
system is defined on a spatial lattice with finite lattice constant also umklapp processes 
are allowed and '0' above can be replaced by any reciprocal lattice vector. 

In the following discussion we will symbolically represent the different interaction 
terms in 1)3111 as follows, 



'^ki,a'^k2,«'^k3,/3^k4,/3 




(33) 



k2 k4 

where spin conservation is understood horizontally. This term represents the process 
where two fermions with momenta k2 and k4 are scattered and thereby obtain new 
momenta ki and ka, respectively. The Hubbard model on a finite 2D L x L lattice 
corresponds to the special case where"*" 

k = {kx, ky), kx^y = {2tx / L) x integer, — vr < k^^y < vr, (34) 

a =t, i, and 

= -2t{coskx + cos ky), v^^ = ^ 5ki-k2+k3-k4,K, (35) 

K 

with t and U > the usual Hubbard constants; the sum is over all reciprocal lattice 
vectors K = {Kx,Ky) with K^y = 27rx integer. Thus in the Hubbard interaction all 

+ To sec that this is identical with the standard definition of the 2D Hubbard model 2 use 
n^ni — : {n-\ + n^)'^ : /2 and perform a lattice Fourier transform. 
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scattering processes occur and have equal weight. The number of distinct 1-particle 
states in this case is 2L^, and the Fermion Fock space is 2^^^-dimensionaL 

We now discuss various simphfied models which can be obtained from a models of 
the type given above. Our method is by truncation, leaving only particular interaction 
terms. This truncation method might seem ad-hoc, and I therefore stress that it only 
should be regarded as a simple way to find the structure of 'interesting' Hamiltonians. 
Systematic derivations of these Hamiltonians should be based on renormalization group 
methods; see e.g. jTHI and references therein. 



3.1. Hartree- Fock models 

As is well-known, not all interaction terms in an interacting fermion model are equally 
important. For example, the terms where the particle momenta are conserved, ki = k2 
and ka = k4, or exchanged, ki = k4 and ka = k2, are known to be of particular 
importance. We will refer to them as as Hartree- and Fock terms, respectively. They 
can be represented symbolically as follows. 



q k q 





+ \/ (36) 



k q k q 

Hartree terms Fock terms 

The simplest example of an exactly solvable model is obtained by truncating the 
interaction in (|31|) and keeping only the Hartree-terms, 

Vn= ^^'^ ■ KcKp ^ (37) 

k,q,a,/3 

where I introduced the fermion number operators 

%,a = 4,aCk,«- (38) 

These operators all commute, and thus the Hartree Hamiltonians H = Hq + are 
sums of commuting terms and can be easily diagonalized. Note that one obtains such 
a model also if one also keeps parts of the Fock-terms, namely those with a = (3: this 
adds to the interaction the terms 

k,q,a 

and all we say below can be easily generalized to Hartree-type models with Hamiltonians 
H = Ho + Vn + V^. 

We naively described the procedure yielding the simplified Hamiltonian H = 
Hq + Vh from the one in Eq. (jHm) -(|M |) as truncation dropping all but certain particular 
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terms in the interaction. This suggests simple formulas of the parameters VI/^k,q in 
terms of the fki,...,k4- However, a systematic derivation should be based on some 
renormalization procedure in which the model parameters are renormalized and all 
but certain interaction terms become irrelevant. Thus this truncation should also be 
accompanied by some change of model parameters. One thus should interpret the 
Hamiltonian H = Hq + Vf^ as a phenomenological model with model parameters to 
be fixed by experimental data. Similar remarks apply to all other truncated models 
discussed below. 

The Hartree Hamiltonian H = Hq + Vhis exactly solvable: as can be checked easily, 
its exact energy eigenstates are given by 

^n = n(i«)"'""^' (40) 
k,o 

and are labeled by the occupation numbers n = {rik.a}, "n-k.a = 0, 1. The corresponding 
energy eigenvalues are 

^n = J2^E^-^'>Ka+ E ^k,qnk,,nq_^. (41) 

k,a k^q,a,/3 

This exact solution makes the analysis of this model much simpler as compared to a 
general interacting fermion model in Eq. (j30|) - (j32|) . but it important to note that this 
solution still is quite far from understanding the physics of this model: to find the ground 
state of the Hamiltonian in ()37|) at fixed particle number one needs to minimize the 
functional in Eq. pT|l over all configurations n such that a ^k,a = N, and even 
though the exact solutions yields the following formula for the partition function, 

Z = J2e-^'-, (42) 

n 

its computation still is a non-trivial task (it is equivalent to solving and Ising-type 
model). It is interesting to note that such an analysis of a Hartree-type model appears 
to be very similar to Landau's Fermi liquid theory [Tij . 

A important feature of the Hartree-type models is that their eigenstates are identical 
with the ones of the corresponding non-interacting model: the interaction only affects 
the energy eigenvalues, but the eigenstates remain Slater states without any fermion 
correlations. A more complicated model taking into account some correlations is 
obtained if one keeps, in addition to Hartree terms, all Fock terms (see fl36|) above). 
This yields an interaction of the following type, 

VhF = E : Wk,q5^k^^q + ^k.qSk " Sq : (43) 
k,q 

with the usual spin- and charge density operators S = {S^, S"^, S^) and n, 

= 5Z4,«(o-)a/3Ck,/3 and ^k = X14,aCk,«, (44) 

a, (3 a 

and cr = (a^, cr^, o"^) the Pauli spin matrices as usual. 



2D interacting fermion models 



14 



The Hamiltonian H = Hq + Vhf is no longer the sum of mutually commuting terms. 
Its eigenfunctions thus are more complicated. To see that it is convenient to write the 
Slater states in Eq. pn|l as follows, 

. . . , a^v) = ■ ■ ■ cl^,aj,^- (45) 
The action of H = Hq + Vuf on these states is easily computed, 

H\ai, . . .,aN) = iSo + ^J{kj,kk)a-j ■ a-j, \ . . ■ .a^) (46) 

V j^k ) 

with Eq = ^j{Ey., — /i) + X]jyfc^kj,kj the Hartree-model eigenvalues and the &j 
Heisenberg spin operators acting on the spin a^-index of the fermions as usual. To 
find the eigenvalues of the Hartree-Fock model one thus needs to diagonalize the spin 
Hamiltonians Ylj^k '^O^j^^k)^j ■ ^k- The latter can be interpreted as ID Heisenberg 
spin systems with long-range interactions. In the case where J(k, q) = J is constant, the 

spin Hamiltonian reduces to J (^f^i o-k^ + const., and its eigenvalues and eigenstates 
can be determined using group theory (i.e. decomposing the direct sum of N spin-1/2 
representations of the Lie algebra of SU(2) in irreps). 

3.2. BCS-type models 

These models underlie the theory of of superconductivity. They can formally be obtained 
from a general interacting fermion model in Eq. ()30 |) - P^ by dropping all but the 
following terms in the interaction. 




-q,i) 

(47) 



(k,T) (-k,i) 
BCS terms 

The resulting Hamiltonian is H = Hq + Vbcs with 

V^Bcs = $^^?k,qA+Aq (48) 

k,q 

where 

A,^=iAa' = ^-q,i^q,T- (49) 

This Hamiltonian is exactly solvable for parity invariant systems where E]^ = -E_k- One 
way to see this is to write Hq = X]k(-^k ~ /^)'S'k (up to an irrelevant additive constraint) 
where 

'5'k = Ck |C^^^ — c_k,j,cLk,|- (50) 
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One can check that the operators A='= and S represent the Lie algebra of SU(2), 

[A+, A^] = [^k, Aj] = ±<5k,qA± (51) 

Thus the BCS Hamiltonian is a quantum SU(2) spin-type model of similar to the 
Hartree-Fock model described above. One can make this relation more specific by 
introducing the following Bogoliubov transformation, 

leaving the canonical anticommutator relations invariant. Then A^ = (S^. ± iS^)/2 
and S'k = with Sk = cj^crcj^ the spin operators associated with these new fermion 
operators. The BCS Hamiltonian can then be written as 

H = Y,iE^.-f^)Sl + Y, 9Kn (^k^i + S'Jl) (53) 

k k,q 

which formally is identical with a xy Heisenberg spin model in a transverse magnetic 
field. This model can therefore be treated in a similar manner as the Hartree-Fock model 
discussed above. The exactly solvable case of constant coupling, g^ ^ = g, underlies the 
BCS theory of superconductivity j3] and was solved exactly in 0. 

3. 3. A different interpretation of the QH model 

Our discussion above suggests the following interpretation the interaction of the QH 
model: we can regard the quantum numbers i and m as the x- and y-components of a 
2D pseudo momentum, i = kx and m = ky. We then can interpret Vqh in Eq. (fT^ as 
mixed scattering terms which are Hartree-like in the x-component and Fock-like in the 
y-component of the momenta or vice versa, i.e. 




{kx, ky) {Qx, (Jy) {kx, ky) {Qx, (Jy) 



mixed terms I mixed terms II 

As discussed, one can exactly solve the generalized QH model H = Hq + Vqh + Vr 
defined in Eqs. (ITnjl . (|T^ and (fTTjl . Interpreting Em + Ei = as dispersion relation 
we can regard this as simplified model obtained by truncating a spin-less variant of a 
model defined in Eqs. (j30|l - (j3T|l . keeping not only Hartree- and Fock terms but also the 
mixed terms just described. One main message of this paper is that it is possible to keep 
these mixed terms and still have an exactly solvable model. 

As is well-known, spatial correlations are important in 2D, and the simplified models 
models including only the Hartree- and Fock terms ()36|) are not adequate. The main 
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suggestion is to use the generalized models in 2D where, in addition to the Hartree-Fock 
terms also the mixed terms ^54\ l are included. 

It is important to note that, in our present interpretation of the QH model, 
interchanging ^ ky amounts to changing ^ Ei and g ^ —g; see Eqs. fll3|) - 
Thus, for spin- less fermions, one can have a non-trivial correlation interaction Vqh 
only if parity invariance is broken. As discussed below, this no longer is the case for 
fermions with spin. 



3.4- Generalized Hartree model 

I now define a parity invariant variant of the QH model including spin. It corresponds 
to a generahzation of a Hartree model Hq + Vn defined in Eqs. fjHIH) and (jTfjl . We 
restrict ourselves to 2D dispersion relation of the following form, E]^ = e^^ + e^y, and 
assume that the Hartree coupling can be written as W^k.q = Vk^q^ +^kyqy +Wkxqy +Wq^ky- 
Introducing the operators 

PkyQy ^^^kxkya^kxQya^ PkyQy ^ ^ ^kr.k,ia^PT.k,ia' (^^) 
kx,a ky,a 

we can write 

Ho = ^ek{pkk + Pkk - p) (56) 

k 

and 

Vh = ^ Vk,q{PkkPqq + PkkPqq) + Wk,q{pkkPqq + PkkPqq)- (57) 
k,q 

I now propose the model H = Hq + + Vmixcd with the following mixed interaction, 

Knixed = ■ ^P'^<lP<lk + PkqPqk ■)■ (58) 

k,q 

This model is obviously parity invariant, and due to the presence of spin the interaction 
T^iixed is non-trivial. It is straightforward to generalize the discussion in 12. 31 and 12. 41 One 
finds that this model also has a dynamical symmetry gl^^ © gl^, and explicit formulas 
for its eigenstates and eigenvalues can be obtained as well. It would be interesting to 
explore the physics of this model in more details. 



3.5. On the 2D Hubbard model 

To illustrate the flexibility of the truncation procedure advocated above, I now discuss 
various simplified models for the 2D Hubbard model. I also show that the exact solution 
of the Hartree-like models for the 2D Hubbard model reproduce mean field theory, as 
claimed in the introduction. 
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3.5.1. Simplified models. For the 2D Hubbard model the Hartree-hke model Hi = 
Ho + Vu + Vf defined in Eqs. dSHl), (EZI) and is given by 

Ho = J2 ^k^k, Vn + VF = ^n^h^ (59) 

k 

with Eif_ in ()35|) and Ua = X]k^k,a- As I show below, this model is equivalent to mean 
field theory (= Hartree-Fock theory) restricted to ferromagnetic (F) states. It misses 
an important physical property of the 2D Hubbard model which, at half filling, has 
anti- ferromagnetic (AF) rather than F order; see e.g. 0. It is interesting to note that 
there is also an exactly solvable model accounting for AF: It it obtained by including 
Hartree-like terms where the momenta are not fully conserved but changed by the AF 
vector Q = (tt, tt), ki = k2 + Q and ks = k4 + Q, and similar Fock-like terms.* These 
yields the following interaction term, 

Vh'^ + ^f'' = ^^^^'^^^ ^^ = E^Ua„- (60) 

k 

The mo del H2 = Hi + + can be solved exactly, and its solution yields Hartree- 
Fock theory for the 2D Hubbard model allowing for F, AF, and charge- density waves. 
The corresponding mixed terms are 

Kiixed = ^ : Pkq^Pqk^ + Pkq^Pqk] '■ (61) 



L2 

k,q 



with pkq and pkq in and similarly V^^^g^. It would be interesting to determine the 

''mixed ' 



phase diagram of the model H3 = H2 + Knixed + K^xed- ^ hope to come back to that in 



future work. 

3.5.2. Relation to mean field theory. I now show how to compute the partition 
functions of the Hartree-like models above, and that this reproduces standard mean 
field theory. For simplicity I restrict myself to the Hartree-like model allowing for AF 
only, 

^ = ^0 + ^(^T + + ^(n^ - (62) 

(we can ignore normal ordering here since = n^, and thus normal ordering only 
amounts to a shift in the chemical potential). It is easy to extend my argument to the 
other models discussed above. We use the following identity (Hubbard-Stratonovich 
transformation) , 

e-^^ = const. / drdse-^(^'(^'+^')/^+^«+^("?-"?)+"-("T+«i)), 



which linearizes the interaction at the cost of introducing two integrations. With that 
we can compute the partition function (we drop the irrelevant constant), 

Z = [ drd.e-^^'[(^'+^')/^-^('-'^)l (63) 
* These are umklapp processes included in the Hubbard interaction in Eq. 13511 . 
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with e ^ the partition function of non-interacting fermions coupled to constant 
external fields r and s. The latter can be easily computed, and one obtains 



the fermion bands in an AF background. Note that is a Riemann sum converging 

to an integral in the thermodynamic limit L —>■ oo. In this limit a saddle point evaluation 
becomes exact. The saddle point equations one thus obtains are identical with Hartree- 
Fock theory for the Hubbard model restricted to AF states; see e.g. 0. 

It is interesting to note that Eq. was obtained by a different method in Ref. |15j . 
In this work we proposed it as a useful alternative formulation of mean field theory for 
the 2D Hubbard model. Superficially this seems equivalent to the 'standard' formulation 
which only uses the saddle point equations resulting from that integral, but it is in fact 
more general since it allows for the possibilities of having degenerate saddle points. We 
found that degenerate saddle points occur in a large part of the parameter regime: if one 
fixes the fermion density close but away from half filling one typically needs to adjust 
the chemical potential /i such that the integral is dominated by two distinct saddle 
points at the same time, one describing AF order, and another describing no order at 
all [T3] . These mixed phases are missed if one only looks at the saddle point equations: 
the mean field phase diagram of the 2D Hubbard model thus obtained is much richer 
than generally believed (see e.g. and references therein). Mixed phases are typical in 
Hubbard-like models. I feel this is a basic, important property of these models which 
deserves to be wider known. Anyway, the present paper gives models for which the 
mixed phases are exact. 

4. Conclusions 

In this paper I advocated the following strategy for treating difficult many-body 
Hamiltonians for interacting fermions: rather than solving such model in some 
approximations (like mean field theory), truncate it so as to obtain a simplified model 
which can be solved exactly. I recalled examples from condensed matter physics a long 
time ago where this strategy was used successfully to solve important problems like 
superconductivity or magnetism. I also discussed a many-body model for 2D quantum 
Hall systems proposed and solved in [T]. I then argued that this model is the simplest 
example in a large class of similar models which should be useful for shedding light on an 
interesting problem in modern condensed matter physics: to understand 2D correlated 
fermion systems. As indicated, the partition function of these generalized models can be 
represented as matrix models. The computations of the latter are interesting problems 
for the future. 




(64) 



with 



£k,± = (^k + ^k+Q)/2 - ± ^s2 + (Ek-^k+Q)V4 



(65) 
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Many of the results reported here can adapted straightforwardly to boson systems. 
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